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The J — > oo double exchange model is formulated in terms of three auxiliary particles. A slow 
true bosonic magnon propagates by admixture with a fast fermionic pseudo-magnon. This process 
involves the absorption of a conduction electron which, for this half-metal, carries only charge degrees 
of freedom. The magnon dispersion becomes much weaker and the relaxation rate increases rapidly 
upon approaching the zone boundary. That the magnons relax for all wave vector values implies 
the existence of a low energy spin continuum. 
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The double exchange model has a long history. It was 
introduced by Zener[l| fifty years ago in order to explain 
the phase diagram of Lai-^Ca^MnOa. With the explo- 
sion of interest in the colossal magnetoresistance effect in 
a number of manganite systems this has again become a 
model of prime importance. The Hamiltonian is, 

H= -Jj2 § i ■ g i - E tnn>(cLc n ,„ + H.C.). (1) 
i crnn' 

Here J reflects the Hund's rule coupling energy and is 
taken to be positive and large, i.e., J > t. In a semi- 
nal paper, de Gennes[^[ |3) showed for near neighbor t nn t 
couplings that in the classical, i.e., large spin 5, limit this 
model has magnon like excitations ex (xt/S)(l — cosg), 
i.e., the dispersion relationship is that of a Heisenberg 
model. For such large values of J the model is a half- 
metal. The Fermi level is in the majority band and Stoner 
excitations have large energies ~ J and are well above the 
much lower energy magnons. 

The experimental situation has been discussed 
recently by Dai et al. Materials with high Curie tem- 
peratures do have Heisenberg like dispersion relation- 
ships but for similar compounds with smaller such tem- 
peratures there is a marked softening of the magnons 
upon approaching the zone boundary. The relaxation 
rate is also anomalous. The magnon width is quite small 
for small q increasing very strongly for q larger than some 
critical value. Since this behavior is not explained by 
double exchange model using the standard theory based 
upon the large S approach nor by the theory of magnon- 
magnon interactions, these authors suggest that the cou- 
pling to phonon modes is important. The theoretical 
understanding has been summarized H by Khaliullin and 
Kilian. They also emphasize that the experimental situa- 
tion cannot be explained by the usual theory and develop 
a large S approach which includes coupling to Jahn- Teller 
phonons and which can explain experiment. 

In this Letter will be presented a novel analytic the- 
ory of the basic TL, Eqn. (jl]), which shows that the two 
experimental effects attributed to phonons do in fact oc- 
cur within the simple double exchange model. This new 
approach is based upon the auxiliary particle method. 



There are three such particles, a 6-boson which reflects 
the true, slow, magnons, the c-fermion which represents 
the charge degrees of freedom of the conduction electrons 
and a t-fermion which corresponds to a spin deviation at 
a site with a charge carrier. 

In fact, the earlier exact diagonalization studies of 
Zang et al., Q and Kaplan and Mahanti [Q for the one 
dimensional model indicate the existence of zone bound- 
ary softening. In particular for S = 1/2 and small con- 
centrations the magnon dispersion deviates strongly from 
(1— cosg). Very recently Kaplan et al., || have performed 
new exact diagonalization studies again for one dimen- 
sion and in the limit J — > oo. Contrary to popular belief 
and mean field theory, they find spin excitations with 
energies much less that J. These they associate with a 
non-Stoner continuum and speculate that this continuum 
comes down to zero energy in the thermodynamic limit. 
In contrast, Golosov || using the large S expansion up 
to fourth order claims that magnon-electron scattering 
does not give rise to magnon damping this implying that 
the magnons do not enter a continuum at low energies. 

Mathematically the present approach is based upon an 
expansion in x the concentration of carriers. For small 

S = | and one dimension our results for the magnon 
dispersion agree well with those of exact diagonalization. 
Useful results are obtained for x as large as 3/16 even 
for S = 1/2. The magnons can exchange energy and 
momentum which the charge excitations and all finite 
q excitations have a non-zero relaxation rate confirming 
the existence a low energy non-Stoner continuum. For 
q < kp : the Fermi wave vector, relaxation is reduced by 
phase space considerations. A rapid onset of relaxation 
will therefore occur around q « kp- Although most of 
the results presented here are for one dimension, these 
principal conclusions are also valid for higher dimensions 
and some such results will be given. The present ap- 
proach does not include orbital degeneracy nor phonons 
and the authors recognize that these are important phys- 
ical ingredients for real materials. The method is easily 
generalized to include such effects and might form the 
basis for a comprehensive theory at least for modest x 
values. 
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Consider a single spin deviation relative to the ferro- 
magnetic ground state for an odd number of electrons. 
In the J — > oo limit, for a given site n, there are only 
four states of interest. There is the state 1*5,0) with 
S z = S and no electron and which is considered to be 
the "vacuum". The state \S — 1, 0) has the spin deviated 
and is mapped to a state d n \) which has a d-boson cre- 
ated on this vacuum. An un-deviated state with an up 
conduction electron is |5, |) and maps to the c-fermion 
state c n \). Lastly, a state which has a total spin S 



but S z = S - : 
t n \) as its map 



is 



2<?|S-1.T) + |S,I) 



The particles are "hard core" 
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with 



total number Q n = n n + + n n < 1, where 



i.e., the 
e.g., 

n n — Ai^n- The total number of conduction electrons 

N = Y^n ( U n + n n) °^ oes not involve d n . 

Subject to Q n < 1 the low energy part of the Hamil- 
tonian maps exactly to, 



1 



2S + 1 



-tl,t n + H.c.) 



2S 



2S- 



l^2tnn>(tld n ,C n + c{,d n t n + H.C.). (2) 



where, e.g., c n = c n (l — n n — n^). When acts on a state 
which satisfies Q n < 1, the result is also consistent with 
this constraint, i.e., this representation of the Hamilto- 
nian does not mix physical with unphysical states. 
The Hamiltonian is written as H! — fiN + I c where 

ic = tj: n [(cic n - ci +1 c n ) + (i/(25 + m^+ik - 

t\i+\t n ) +H.c] comprises the operators which enforce the 
constraint. Then the Fourier transform of, 

H'-nN = 



25 



26* 1 

- — V eMt J t Cff + cUU r , _ + H.cX (3) 

k,q 



where = — . tije lk '( rj rj \ The most important part 
of I c is, 

iV^-fct c r dt A^+^.c, (4) 

^2 <*■ \ k+q k k'-q k' J ' w 



A/" 2 

since this leads to the renormalization of the conduction 
electron wave function near the deviated spin. In the 
magnon dispersion relationship, the leading term 0(x,i) 
is an unimportant constant. The 0(x,t 2 ) term produces 
renormalization and relaxation and will be discussed be- 
low. 

In the zeroth order approximation, for small x, the 
conduction electrons occupy the region near k = at the 
bottom of a band of total width D. The t-particle band 




FIG. 1: For one dimension, the conduction electron, or c- 
particle, band has width D = At while the pseudo- magnon, 
i.e., t-fermion, band is a factor of 2S + 1 narrower 




FIG. 2: (a) A true, d-particle, magnon decays into a pseudo- 
magnon, i.e., a t-particle, upon absorbing a conduction elec- 
tron c-particle. The momentum is conserved and each vertex 
is y/2S/(2S +T)(e r /JV). This leads to Eqn. %. (b) The 
equivalent for the ferromagnetic state. The bo -boson, which 
reflects the state 15,0), decays to a /o-fermion, which corre- 
sponds to \S, T), upon absorbing a conduction electron. This 
is reflected by Eqn. (g). 



has a smaller width [D/(2S + 1)] and for small x this 
band will not have a real occupation, see Fig. 1. The 
deviated spin, associated with the d-particle, has no zero 
order dispersion. This becomes dispersive by admixture 
into the t-particle. There are matrix elements, y/2S, of 
the total spin which connect the vacuum with the boson 
d\, i.e., this is a true magnon. There are also matrix 
1, which connect c k with t\ +q , so the 



elements, \/2S~- 

finite wave- vector total spin operator, 



254 + V2STT 5>g +f c E + H.c). (5) 



For small x the dispersion in energy and momentum of 
the conduction electron, c-particles, is relatively small 
and so, for the second term, the t-fermion dominates the 
q dispersion. Since it has a band width ~ D/S the t- 
fermion can be identified as the fast pseudo-magnon with 
a large gap ~ D/2. 

The leading approximation for the dispersion of the 
slow magnon, corresponds to the Feynman diagram 
shown in Fig. 2a, 



2S 



u d > = 



2S- 



1 N 2 — ^ k 



2S+l € k+q 



(6) 



where is the usual conduction electron thermal oc- 
cupation factor. Since it involves a single sum over the 
occupied conduction electron states uo^^ x which shows 
that the <i-true-magnon moves slowly for small x. 



3 



It is important to confirm that the present method 
does not give a false spin gap. In order to compare like 
with like, the problem with no reversed spins is formu- 
lated using a method in which a single site is treated 
specially so that the same approximation as used above 
for the deviated <i-particle site, can be made for this spe- 
cial site. There are two, rather than a single operator 
for this "impurity" site, taken to be at the origin,. The 
absence of a particle corresponds to an unphysical vac- 
uum |)o. The maximum spin state without a conduction 
electron, \S, 0), is mapped to Do while \S, |), when a 
conduction electron is present, is reflected by /q|)o- The 
manifold with no reversed spins then corresponds to, 

nn' n^O 

+ ^ton (f^b (c n ^c N _ n+1 )^H.c. -/i/oVo- (7) 

With this unusual representation of the ferromagnetic 
ground state it is possible to make the same approxima- 
tions as with a single deviation and a self-energy, Fig. 2b, 
for the 6-particle, the equivalent to the d-particle. This 
gives for the energy, 



1 2 n k = 1 

k k 



(8) 



which, surprisingly, is the exact answer. Clearly Eqn. (g) 
agrees with Eqn. (^) when q = which indeed confirms 
that there is no gap in the magnon spectrum. 

It should also be that uo^ agrees with the (corrected) 
result of de Gennes for large enough S. The large S 
expansion for the one dimensional near neighbor model, 
i.e., with eg — > —2tcosk gives, 



1 

2S 



N ^ 



cos fcru 



(1-cosg) (9) 



which agrees with FurukawaQ. 

It is also possible to obtain an analytic result for small 
enough x independent of S. Assuming that the conduc- 
tion electron Fermi surface contracts to a point, gives, 



U) 



4xtS 



25 + 1 1 



(10) 



25+1 



For small values of 5, this expression gives a dispersion 
which becomes rather flat near the zone edges, see below. 

It is the slow <i-magnon which is to be identified with 
the low energy magnon branch found in the numerical 
work [[I [t[ ||]. In Fig. || are compared the numerical 
results || for S = 1/2 with the present Eqn. (||). For 
x = 1/16 and 3/16 the flattening of the dispersion near 
the zone edge seen in the numerical data is well repro- 
duced although the difference between the two different 




FIG. 3: The right box corresponds to the present theory 
for concentrations x = 1/16 and x = 3/16 for spin S = 1/2. 
The top curve corresponds to A(l — cosq). The scales where 
chosen so that all curves coincide at q = 0.1257T. The small x 
result Eqn. (HQ) cannot be distinguished from the result for 
x = 1/16. The left box shows the numerical results || with 
n electrons and 16 sites. Again the scales cause the points at 
q — 0.1257T to coincide with each other and the A(l — cosq) 
solid line. 



concentrations is not as great as in the numerical data. 
For 5=1/2 and x = 5/16 the Fermi surface lies in the 
t-particle band and the assumptions used for Eqn. @ 
are not valid. The result Eqn. and the small x and 
the large S approximations are compared, for S = 1/2, 
in Fig. [|. The small x approximation compares well with 
x = 1/16 even for S = 1/2. In contrast the agreement 
with the large 5, simple cos A:, result is not at all good. 
The deviations from a simple cosine are therefore well 
reproduced by Eqn. (|io| ) and are clearly most extreme 
for S = 1/2. 

As the small x approximation illustrates, the fatten- 
ing near the zone edges occurs because, for small 5, the 
dispersion of the virtual t-particle is comparable with its 
energy. 

Since the effective exchange increases with increasing x 
so will the Curie temperature, T c . Larger x values, Fig. ||, 
have a dispersion which approaches that of a Heisenberg 
model. It follows for larger x and T c values, the small q 
magnon stiffness will accord with T c . It is a commonplace 
observation that the zone edge softening which occurs 
for smaller x values will cause T c to fall short of that 
predicted by this initial stiffness. This appears to be the 
case experimentally Q. 

It is evident from the Hamiltonian, in the form 
Eqn. (||), that there is a continuum of low energy longitu- 
dinal, i.e., charge excitations of the conduction electron 
c-particles. The question is whether the spin excitations 
couple to the continuum of charge excitations. This is the 
case and such a coupling shows up as a finite relaxation 
rate for the d-magnons. In one dimension the leading or- 
der relaxation process for this slow magnon has an onset 
for relaxation which occurs only at finite momentum in- 
dicating that this low energy branch enters a continuum 
for finite momentum. This relaxation corresponds to the 
Feynman diagrams shown in Fig. Corresponding to 
Fig. ||a is a second order process which involves J c , while 
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FIG. 4: Curve "3" is Eqn. ([]), with x — 0.15, for u q and 
is almost indistinguishable from curve "2" which reflects the 
approximation - — — 2xt. Curve "4" is the large 5 

approximation 4a:t[S'/(5 r + 1) 2 ](1 — cosg) while curve "1" is 
similar but made to fit the q 2 dependence of "3" for small q. 

Fig. ||b is of fourth order and uses the vector t-fermion in 
the extra intermediate states. There is also interference 
between the two processes. The net relaxation rate is of 
the form: 

t' 2 

p',q 

x5(uo q + 2t cos k — (uj p + 2t cos(q + k — p)) , (11) 

where t f « t + [25/(25 + l)](2t 2 /|/i|). 

A particle-hole charge excitation of the conduction 
electron is created and takes away energy and momentum 
from the magnon. Since only the initial state electron 
lies below the Fermi surface, such a process is propor- 
tional to x. The only energy scale is t (e.g.,|/i| ~ t) but 
relaxation is strongly suppressed due to phase space con- 
siderations. In order to understand the trends, consider 
small x and larger values of 5. This implies that the en- 
ergy dispersion of the magnon is negligible as compared 
to that of the conduction electrons. It is not possible to 
conserve both energy and momentum for small values of 
the magnon q since the energy of the conduction elec- 
tron charge excitations with the same momentum is too 
large. Small energy but large q particle-hole excitations 
imply that the electrons scatter from ^kp — > ±kp and 
have a momentum transfer of ~ 2k p. In order to mini- 
mize the initial magnon momentum it must scatter from 
^kp — » ±fcp, i.e., the onset of magnon relaxation will 
occur for q w ±kp a momentum which is proportional to 
x for one dimension. The finite magnon dispersion im- 
plied by smaller values of 5 and larger x will reduce this 
threshold value of q from this limiting value of ±kp. 

However higher order processes have a final state 
which still comprises a single magnon but now reflects 
a pair of conduction electron particle-hole excitations. 
Even with an arbitrary small initial magnon momentum 
it is always possible to find a pair of particle hole exci- 



tations which, e.g., go from close to one Fermi surface to 
the other (±kp — > =F^f) but in opposite senses and which 
have any small energy and moment. Thus as x increases 
the pseudo-gap to momentum « ±kp progressively fills. 




FIG. 5: (a) A momentum q magnon relaxes via I c to momen- 
tum p transferring energy and momentum to a particle-hole 
conduction electron pair, (b) The same process but using the 
t-fermion as intermediate particle. 

Much of the above discussion generalizes to two and 
three dimensions. In particular the small x approxima- 
tion, 

_ 2xe- k=0 S 1 1 
9 25+1 Nl-^e./e^) 1 ) 

exhibits softening at the zone edges and relaxation will 
be phase space limited for q < kp even though there 
is no longer a threshold for the second order process. 
Evidently the orbital degeneracy of the real materials is 
of great importance. This will be dealt with elsewhere. 
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